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Abstract 

In this paper we point out that the cosmological baryon asymmetry in our 
universe is generated almost independently of the reheating temperature Tr in 
Affleck-Dine leptogenesis and it is determined mainly by the mass of the lightest 
neutrino, niy^, in a wide range of the reheating temperature Tr ~ 10^-10^^ GeV. 
The present baryon asymmetry predicts the rriy^ in a narrow region, ~ (0.3- 
1) X 10~^ eV. Such a small mass of the lightest neutrino leads to a high pre- 
dictability on the mass parameter ruy^y^ contributing to the neutrinoless double 
beta decay. We also propose an explicit model in which such an ultralight neutrino 
can be naturally obtained. 



1 Introduction 



The origin of baryon (matter-antimatter) asymmetry in the present universe is one of 
the fundamental problems in particle physics as well as in cosmology. Recently, lepto- 
genesis becomes very attractive among various baryogeneis scenarios, since there are 
now convincing evidences of neutrino oscillations, especially the atmospheric neutrino 
oscillation reported by the Superkamiokande collaboration p|. The small but nonzero 
masses of neutrinos suggested from neutrino-oscillation experiments strongly indicate the 
existence of lepton-number violation, which is a crucial ingredient of the leptogenesis. 
It is extremely interesting in the leptogenesis scenario that the baryon asymmetry of 
the present universe is closely related to observable low-energy physics, namely, neutrino 
masses and mixings. 

Among various mechanisms @, @, H, H for leptogenesis proposed so far, the lep- 
togenesis of Affleck-Dine (AD) mechanism is naturally expected to work once one 
introduces supersymmetry (SUSY) in the standard model together with the nonzero 
neutrino mass. Recently, a detailed analysis on this mechanism (AD leptogenesis |^) 
has been performed [§, which has shown that the dynamics of the flat direction fleld is 



drastically changed by thermal effects, as pointed out in Refs. [£], |T0|. Actually, it has 
been shown [§] that the resultant baryon asymmetry is indeed suppressed for relatively 
high reheating temperatures Tr, and an ultralight neutrino is required in order to obtain 
the sufflcient baryon asymmetry hb/s in the present universe. (Here hb and s are the 
baryon number and entropy density in the present universe, respectively.) 

In this paper, we perform a reanalysis on the AD leptogenesis, and emphasize that the 
present baryon asymmetry ub/ sis determined mainly by the mass of the lightest neutrino 
uiy^ for Tr ~ 10^-10^^ GeV, and its dependence on the reheating temperature Tr is rather 



mild. Here, we include an additional thermal effect observed recently in Ref . [TT[ , which 
makes the dependence of the resultant baryon asymmetry on the reheating temperature 
even milder. Notice that in many baryogenesis scenarios the obtained baryon asymmetry 
depends crucially on the reheating temperature Tr of the inflation. Thus, it is very 
attractive that the baryon asymmetry in our universe is predicted almost independently 
of the reheating temperature Tr in the AD leptogenesis. 

Furthermore, this reheating-temperature independence of the AD leptogenesis means 
that the mass of the lightest neutrino iriy^ can be determined from the present baryon 
asymmetry ub/ s. Actually, we show that the observed baryon asymmetry ub/s c::^ (0.4- 
1) X 10~^° [|12| predicts m^^ ~ (0.1-3) x 10~^ eV. Thus, neutrinos can not be degenerate 
in mass, and the masses of the two heavier neutrinos are also determined from the 
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neutrino-oscillation experiments for the atmospheric and the solar neutrinos, that is, 
m^3 ~ JSrr^ ~ (3-8) x lO'^ eV and m^^ ~ J^r^i ~ IQ-^-lQ-^ eV ]13 . 



Although it is hard to confirm the mass of such an ultralight neutrino Vi, it can be 
tested indirectly. A crucial observation here is that such a small mass of the lightest 
neutrino together with the masses and mixings of the neutrinos obtained from neutrino- 
oscillation experiments suggests a high predictability on the rate of neutrinoless double 
beta {Qv(3(3) decay. We show that the Ve^i^e component of the neutrino mass matrix, 
^ueUe, contributing to the OuPP decay, is determined with high accuracy, depending 
on the solution to the solar neutrino deficits and the e-3 component of the neutrino 
mixing matrix, Ues- For the case of large angle MSW solution, which is favored from 
the recent Superkamiokande experiments |]l3l, a sizable m^^^^, say, \mi,^^J ~ 10^^-10^^ 
eV is predicted. It is very encouraging that such a large \mu^uj is accessible at future 
experiments for the Oi//3/3 decay such as GENIUS fl^. On the other hand, we find 
that the obtained \mu^uj depends highly on Ues for the case of small angle MSW and 
LOW solution. We also stress that this predictability on the Oz//3/3 decay is a generic 
consequence of the mass hierarchy m^-^ <^ m^^j 3 . 

Finally, we propose an explicit model based on a Froggatt-Nielsen (FN) mecha- 
nism in which such an ultralight neutrino is naturally predicted. Here, we impose 
a discrete Zq group as the FN symmetry. 



2 Affleck-Dine leptogenesis 

The Affleck-Dine (AD) flat direction scalar field for leptogenesis is 

Hu = Li = , (1) 

and along this direction we have effective dimension-five operators in superpotential; 

W = ^ (L,hS = (LiHS , (2) 

where = sin/3 x 174 GeV and tan/? = {Hy) / {Hd}- Here Hu (Hd) and Li are 

Higgs field which couple to up (down) type quarks and SU(2)j^-doublet lepton fields, 
respectively. (They denote chiral superfields or their scalar components.) Hereafter, we 
take sin/5 ~ 1. Note that the scale Mj are related to the neutrino masses m^,. through 
the seesaw formula rriu^ = (Hu)"^ /Mi [|T6|. Here we have taken a basis where the mass 
matrix for neutrinos Ui is diagonal, and for simplicity we will suppress the family index 
i (= 1,2,3) unless we denote it explicitly. 
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Before discussing the detailed dynamics of the field, we first roughly describe the 
evolution of and note the relevant epoch for the AD leptogenesis. After the end of 
inflation, the infiaton x starts a coherent oscillation. (At this stage the energy density of 
the universe is dominated by the oscillating infiaton x ^"^^ the Hubble parameter H of 
the expanding universe decreases with cosmic time t as H = (2/3)^"^ pl.) After that. 



when the Hubble parameter H becomes comparable to the decay rate of the infiaton, 
~ T^/M^ (M* = 2.4 X 10^^ GeV is the reduced Planck mass), the energy density 
of the radiation starts to dominate the universe. The evolution of the field is as 
follows. During the inflation, the field takes a large value determined from the effective 
potential discussed below. After the end of inflation, the value of gradually decreases 
as the Hubble parameter H decreases and then, at some time, the starts its coherent 
oscillation. As we will see, the net lepton number is fixed when the flat direction field 
starts the coherent oscillation. 

Let us discuss the dynamics of the field. The method of the following analysis in 
this section is based on Ref. [§. First, we show the total effective potential for the flat 
direction field relevant to the leptogenesis. In addition to the usual F-term potential 
and SUSY-breaking terms, there are additional SUSY-breaking terms caused by the 
nonzero energy density of infiaton, which depend on the Hubble parameter H, and also 
there are thermal potential terms depending on the temperature T. It turns out that 
the total potential is given by the following form;Q 

V fk\<t>\<T 




We explain each terms by turns. First of all, the F-term potential directly comes from 
the superpotential Eq. (g); 

v.=^i^r (4) 

Next, the SUSY-breaking potential for is given by 

SVo = mj|0p + ^ (a^0^ + H.c.) , (5) 



^Thc thermal potential proportional to T"* log in Eq. which has been recently found out 

in Ref. Illll, was not considered in Ref. H. 
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where and m3/2am are SUSY-breaking mass parameters. They are expected to be 
~ "^3/2 ~ 1 TeV and \am\ ~ l.Q 

During the inflation and during the oscillation period of inflaton x after the end of 
inflation, the energy density of the universe is dominated by the inflaton x- Thus, there 
appears additional SUSY-breaking effects caused by the nonzero energy density of x 0; 

5Mnf = -chH^\<P\^ + ^ (aH0' + H.c.) , (6) 

where ch and a^j are real and complex constants, respectively. Hereafter, we assume 
ch — {cihI — 1; since we find that the conclusions in the present paper do not depend 
much on these parameters unless ch ^ 

The rests in Eq. (^) correspond to the thermal effects which we discuss now. Although 
the energy density is dominated by the inflaton x during the inflaton-oscillation period, 
there is a dilute plasma consisting of the decay products of the inflaton x even in this 
period. The temperature of this dilute plasma is given by [|T^ 

T={T^M,Hy^\ (7) 

This dilute plasma has crucial effects on the dynamics of the field. First, the fields ipk 
which couple to are produced by the inflaton decay and/or by thermal scatterings if 
their effective masses are less than the temperature; 

M0|<T, (8) 

and hence the flat direction field receives a thermal mass of order ~ f^T. Here, fk 
denote Yukawa or gauge coupling constants of ipk to 0. The induced thermal mass term 
is given by [p!0| ; 



E c,f,T'\<p\\ (9) 

where are real positive constants of order unity. (Details for and couplings fk 
relevant to the flat direction can be seen in Ref. ||^.) 

Moreover, it has been recently pointed out [jTT| that there is another thermal effect. In 
the 0/v^ = L = Hu flat direction in Eq. (|^), the SU(3)c gauge group is not broken, and 
gluons, gluinos and down-type (s) quarks remain massless. These fields generate a free 

^In this paper we assume the gravity-meditation model of SUSY-breaking. 
^See, for example, Ref. ]l7|. 
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Figure 1: A schematic behavior of the SU(3)c gauge couphng a^. The dashed hne 
represents the running couphng when the (f) field does not have a vacuum-expectation 
value. 

energy which is proportional to gsiT^T^ at two loop order,^ where Qs is the coupling 
constant of the SU(3)c gauge field. On the other hand, the evolution of the running 
coupling Qs is given by; 

_A_,(„._|L(_3C.. E^^T,.,), (10) 

where C2 = 3 and T{Ri) = 1/2 for fundamental representations. Notice that the evolu- 
tion changes when the scale fi passes through an effective mass of a field, fi\4>\, as shown 
in Fig. Ij. In our case denote Yukawa couplings of up-type quarks. Thus, the coupling 
constant of SU(3)c depends on |0| as follows; 

^^(^)Ui>T = ^.(T)i,.o + ^^^E^m^ , (11) 

where i/u are Yukawa coupling constants of up-type quarks and Mq is the ultraviolet 
scale where Qs is fixed. Then, there is an additional potential through the modification 

"^In this flat direction, the down type (s)quarks also generate a free energy of order 0(j/^), where 
yi, is the Yukawa couphng of the bottom quark. This gives an analogous free energy, which does not, 
however, give a dominant effect as long as 1- 
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of the gauge coupling constant in Eq. ([TI|) ; 



5Vt = a,a]T^\ ^ ) log ( ^ ) , (12) 

\yu\'t>\>T 



where ag is a constant which is a bit larger than unity |TT]| and as = g1/{^Ti). Hereafter, 
we take the factor X](2/3)T(i?) to be unity since it does not change the result much.0 

Now we obtain the total effective potential Eq. (|^), i.e., V^otai = Vp + 5Vq + (5V^nf + 
SVi"^ + ^Vj'^. The evolution of (p is described by the equation of motion with this Kotai 
as 

+ 3if0 + ^ = O, (13) 

where the dot denotes a derivative with time. 

During the inflation, there is no plasma and hence no thermal effects. In this stage 
the Hubble parameter H takes a constant value which is much larger than the soft SUSY- 
breaking mass m^. Thus, the potential is dominated by the Hubble- induced terms (5l^nf 
in Eq. (^ and term in Eq. (1), and hence the flat direction (f) rolls down to the 
minimum of the potential. 



101 ~ ^MH 

(A\ -arg(a/^) + (2n + l)7r 
arg(0) ~ , n = 0-3. (14) 

Note that we have assumed the Hubble-induced mass squared is negative {ch — +1 > 0). 

After the inflation ends, the infiaton x starts to oscillate and its decay produces a 
dilute plasma. However, the potential is still dominated by Hubble-induced terms and 
101^ term at the first stage of the oscillation. Thus, the fiat direction field is trapped 
for a while in the above minimum given in Eq. (0) . 

Then, as the Hubble parameter decreases, the negative Hubble-induced mass term is 
eventually exceeded by another term in the potential; 

fk\<t>\<T 

As we shall see below, it is this time when the oscillation of starts. Let us denote 
the Hubble parameter at this time by -ffosc- Using the relations |0| ~ \/MH and T = 



^At least the top Yukawa coupling ft always satisfies /t|0| > T before the oscillation of 0. Thus, the 
resultant baryon asymmetry changes by only a factor of at most. See Eq. dl9). 
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(T^M^H)^^^ , it can be calculated from Eq. (|I5]) as; 



Hr,^r ^ max 



(16) 



where Hi comes from the effect of the coupling /j and is given by 



if.^min|-l^,(c?//M.Tl)^/^}. (17) 

The evolution of the after H ~ i^osc depends on which term in Eq. (|^) dominates 
the effective potential. There are basically three cases; the potential is dominated by 
(i) m^|0p term, (ii) T^|0p term, or (iii) T^log(|0p) term. First, if the potential is 
dominated by the m^l^p term, the equation of motion Eq. ([I3|) is given by 

4> + 3H^ + ml(j) = . (18) 

It is clear that the field oscillates around the origin (0 = 0) and the amplitude of the 
oscillation dumps as |0| oc oc t^^. Second, when the potential is dominated by the 
thermal mass term Cfc/|T^|0p, the field oscillates around = and the amplitude 
dumps as |0| oc H"^^^ oc 0. The third case is given when the a^a^T^ log(|0p/T^) 
term dominates the potential. If we neglect the time dependence of T^, the damping rate 
of the oscillation amplitude due to such a flat potential, V ~ log(|0p), can be estimated 
by using the virial theorem, and it is given by |0| oc H^ oc [|l^]. In the actual 
case, however, the potential itself gradually decreases with time as T"^ oc t~^. We have 
numerically checked that the amplitude dumps as |0| oc if" oc with a ~ 1.5. Notice 
that, in all above cases, the dumping rate is faster than the rate before the beginning of 
the 0's oscillation (|0| oc H^/"^ oc t'^/"^). 

Finally, we derive the resultant lepton asymmetry generated by the AD leptogenesis 
mechanism. Since the field carries lepton charge, its number density is related to the 
lepton number density ul as 



nz. = -2 (0*0-0*0) . (19) 



From Eqs. (H), ([TB|) and (11^, the evolution of is given by 



riL + 3HnL = ^Im (a^0^) + ^Im (0^0") • (20) 

A nontrivial motion of along the phase direction can generate a net lepton asymmetry^ . 
Although the field almost traces one of the valleys in Eq. (|1^, the phase of is slightly 
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kicked by the relative phase between and a/^ in the total potential in Eq. (|^) dur- 
ing its rolling towards the origin. Thus, the first term in Eq. (|20|) plays a role of the 
source of the lepton asymmetry.^y integrating Eq. (pO]), we obtain the resultant lepton 
asymmetry at the time t, 



R^ul] (t) ^ dt'R'lm (a^0^) , (21) 

where R denotes the scale factor of the expanding universe, which scales as R^ oc H~'^ oc 
in the universe dominated by the oscillation energy of the inflaton x- We can see 
that the total lepton number increases with time as R^ul oc t until the oscillation of 
starts {H > i^osc), since 0^ oc and hence ~ const in this stage. On the 

other hand, after the (p starts its oscillation, the production of lepton number is strongly 
suppressed. This is because Im (a^^^) changes its sign rapidly due to the oscillation of 
(j), and also because the amplitude of 0's oscillation damped as R^cp'^ ~ t""^ with n > 1. 
[See discussion below Eq. ([18|).] Therefore, as mentioned in the beginning of this section, 
the net lepton asymmetry is fixed when the oscillation of (p starts. The generated lepton 
number at this epoch is given approximately by 



= (am(p^) t = lm3/2MHosJeS , (22) 



2M 



where S^s — sin(4arg0 + arga^) represents an effective CP violating phase. From 
Eq. (P^, the lepton-to-entropy ratio is estimated asQ 



riL MT, 



R 



C^)S.. (23) 



s 12M2 \Ho 

when the reheating process of inflation completes. This lepton asymmetry is partially 



converted |1| into the baryon asymmetry through the "sphaleron" effects |T9[ , since it is 
produced before the electroweak phase transition at T ~ 10^ GeV. The present baryon 
asymmetry is given by pO|l P| 



s 23 s ^ ^ 

^The contribution to the lepton asymmetry from the second term in Eq. ( |20f) is always less than or 
comparable to that from the first term, since Im(a/f0'*) is suppressed. See Eq. (|l4|). 

^In deriving Eq. ( ^ ) we have assumed that the lepton number is fixed before the reheating process 
of the inflation completes, namely, Hosc^T^/M^, which is satisfied as long as T^^IO^^ GeV for 
< 10"^ eV. 

^In the present analysis we neglect the relative sign between the produced lepton and baryon 
asymmetries. 
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Figure 2: The contour plot of the baryon asymmetries ub/s in the my^-T^ plane. The 
lines represent the contour plots for hb/s = 10~^, 10~^°, 10~^^, and 10~^^ from the 
left to the right. The short-dashed lines represent the baryon asymmetry when one 
neglects the thermal effects. The long-dashed lines represent the ones including only 
the thermal mass oc T^|0p. The solid lines represent the baryon asymmetry including 
both thermal mass and log(|0p) terms. The shaded region corresponds to the present 
baryon asymmetry, n^/s ~ (0.4-1) x 10~^°. We have taken = 1 in this figure. (See 
discussion in the text.) 



Thus, after all, the present baryon asymmetry is give by 

hb _ 2 MTr 



s 69 M2 \H„ 



(25) 



We see that the produced baryon asymmetry becomes larger as the scale M increases, 
i.e., as the neutrino mass rrii, decreases. Therefore, the relevant flat direction for the AD 
leptogenesis corresponds to the first family field, i.e., 4>/V2 = Li = Hu- 

Fig. 1^ shows the contour plot of the produced baryon asymmetry in the tti^-Tb, 
plane. (Here we have used the relation m^, = (Hu)"^ /M.) As shown in the figure, the 
present baryon asymmetry ub/s is determined almost independently of the reheating 
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temperature for a wide range of 10^ ^ ^ 10^^ GeV. In particular, for a relatively high 
reheating temperature 10^ ^ Tr ^ 10^^ GeV, the baryon asymmetry derived from the 
Eqs. (Jl^) and (p5D is given by the following simple from; 

The reason why it is independent of the reheating temperature Tr is that the oscillation 
time Hose is determined by the thermal potential T'* log(|</)p) in the higher temperature 
regime, and Tr dependence is canceled out in Eq. (^5]). Even in the lower reheating 
temperature region 10^ ^ Tr ^ 10^ GeV, where ifosc is determined by the thermal-mass 
term potential T^|0p, Tr dependence is still mild, i.e., ub/s oc T^^^. The reheating- 
temperature independence discussed here is a very attractive and remarkable feature of 
the present mechanism since the produced baryon asymmetry crucially depends on Tr 
in many other baryogenesis scenarios. 

In Fig. I], we have taken 6es = 1- It is expected that 5efr ~ C'(l), say, 6cs — 0.1- 
1, unless there is an unnatural cancellation between arg(am) and arg(ai^). Thus, the 
present baryon asymmetry in our universe n^/s ~ (0.4-1) x 10^^° suggests an ultralight 
neutrino of a mass m^, ~ (0.1-3) x 10"^ eV for Tr ~ 10^-10^^ Q^y ^nd 5eff ^ 0.1-1. We 
consider the region of 10^ < Tr < 10^2 throughout this paper, since it is the case for 
a large class of the inflation models proposed in the framework of supergravity |^T], p^ . 
In Sec. ^ we will propose a model in which such an ultralight neutrino can be naturally 
obtained. 



3 Prediction on the rate of neutrinoless double beta 
decay 

The neutrinoless double beta (OuPP) decay, if observed, is the strongest evidence for the 
lepton number violation. The crucial parameter to determine the OujSlS decay rate is 
Irriiy^i^J = \ J2iU^i'^u,\, where Uai is the mixing matrix which diagonalize the neutrino 
mass matrix.]^ If the mass of the lightest neutrino is actually so small, m^-^ ~ 10~^ eV, as 
discussed in the previous sections, the contribution from m,^^ to nii,^,^^ can be completely 
neglected and hence the parameter \miy^iyj is written in terms of masses and mixings of 
two other neutrinos as 

hvJ_= I Ue2mu2 + U^3^U3 I • (27) 

^As for general studies of the \m^^^J using the neutrino masses and mixings, see, for example, 
Rcf. and references therein. 
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Therefore, the \m^^^J is determined by t/es and the parameters of atmospheric and solar 
neutrino oscillations; Smlf-^ ~ m^^, Sm'^^i ~ m^^ tan^ 6^801 = |?7e2/^eiP- Namely, it 
is given by 

(1 - |[/e3|') Sm^e,ol\/S^l+ |t/e3|V"v^ 



atm 



(28) 



where a denotes the relative phase between the two terms. 

We calculate the predicted value of \m^^^J for the large angle MSW, the small angle 
MSW and the LOW solutions, taking the parameters allowed for atmospheric and solar 
neutrino oscillations that are shown in Fig. ^ and Fig § [24]. In the case of large 



angle MSW solution, \m^^^J is sensitive mainly to the parameter of the solar neutrino 
oscillation, sin^ 9so\\j5ml^y The predicted value of \Tny^i,J\ for the large angle MSW 
solution is shown in Fig. ^. Here, we have required |f/e3| < 0.15 from the CHOOZ 



experiment For a comparison, we also show the possible values of \Tny^^J\ when we 
allow m^^ to be relatively large as mj,^ < {1/ \/2)my^. We see from Fig. | that the \my^v^ \ 
is predicted in a narrow range. It is very encouraging that the predicted \m^^^^ \ can be 
large enough to be accessible at future experiments such as GENIUS |1^. Furthermore, 



if the f/es becomes more constrained by future experiments |^, \my^^J\ is predicted in a 
much narrower range as shown in Fig. where we have required |?7e3| < 0.10. 

On the other hand, the \m,^^,^J\ is sensitive to |f/e3| in the case of the small angle MSW 
and the LOW solutions. The results are shown in Fig. |^ and ||. Because \Uez\ is highly 
constrained by the CHOOZ experiment, the predicted value of \my^yj\ is so small. Even 
in these cases the contribution from m^,^ can not enhance \my^y^ \ because it is too small. 

In all cases, the \my^yj\ is predicted with high accuracy, depending on the solar and 
atmospheric neutrino oscillation parameters and Uez- Therefore, the presence of such an 
ultralight neutrino indicated from the present baryon asymmetry can be tested at near 
future experiments. However, notice that the results shown in this section is a generic 
consequence of the mass hierarchy m^^ <^ "^1^2,3- Thus, we consider that the Oi//5/3 decay 
provides only a consistency test for our hypothesis m^j^ ~ 10~^ eV. 



4 A model for the ultralight neutrino 

In Sec. we have shown that the baryon asymmetry in the present universe predicts the 
mass for the lightest neutrino in a narrow region, iriy^ ~ (0.1-3) x 10~^ eV. Together 
with neutrino masses required to solve the solar and atmospheric neutrino anomalies, 
this suggests a very large mass hierarchy between the lightest and the heavier two neu- 
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trinos. In this section we show an exphcit model based on a Froggatt-Nielsen (FN) 
mechanism in which such an large mass hierarchy is naturally obtained. 



We adopt a discrete as the FN symmetry instead of a continuous U(1)fn- We see 
that the discrete symmetry is crucial to produce the required large mass hierarchy in 
the neutrino sector. Before describing our model, we briefly review the FN model [|T5[] , 
which explains the observed hierarchies in quark and lepton mass matrices. This model 
is based on a U(1)fn symmetry that is broken by the vacuum-expectation value of 
($) 7^ 0. Here $ is a gauge singlet FN field carrying the FN charge = —1- Then, all 
Yukawa couplings arise from nonrenormalizable interactions of $ and are given by the 
following form; 

= h,, e«'+«^ ^^^JH^id) , (29) 

where hij are 0{1) coupling constants, Qi are the FN charges of various chiral superfields 
and e = ($) /M*. The observed mass hierarchies for quarks and charged leptons are 
well explained by taking suitable FN charges for them. For instance, we assign FN 
charges (a, a, a + 1) for lepton doublets Li, while giving charges (0, 1, 2) to the right- 



handed charged leptons Ei, with e ~ 0.05-0.1 Here, we take a = or 1 



Charges for the quarks are determined if one assumes the SU(5) grand unified theory 



The mass matrix for the neutrinos in this model is determined by the FN charges for 



the lepton doublets Lj |23]. To see this we first discuss the mass matrix for the heavy 



right-handed neutrinos iVj, which is given by; 

M^,, = (7., e^'+'^^Mo, (30) 

where Mq represents some right-handed neutrino mass scale and Qij are coupling con- 
stants of 0{1) like hij. Hereafter, we will take a basis where the mass matrix for the 
charged leptons is diagonal.^" The charges for the lepton doublets Li and right-handed 
neutrinos Ni are listed in Table. |I|. Then, the neutrino Dirac mass matrix mu and the 
right-handed neutrino mass matrix Mpt are given by the following forms; 

mo = (Hu) 










( hn 


hi2 


h 







I] 




h22 


h 












h2 


h 




""^"One might wonder if the mixing matrix from the charged lepton sector would change the discussion 
above, since the mass matrix for the charged leptons has off-diagonal elements in the FN mechanism. 
However, the correction from this effect yields higher order terms in e, and hence we can safely neglect 
it. 



12 



Li N3 N2 Ni 



a + 1 



d 



Table 1: The FN charges of lepton doublets and right-handed neutrinos, a = or 1. We 
assume, for simplicity, < b < c < d. 





/ 









( 9u 


912 


913 


Mr = Mo 











912 


922 


923 




v 










\ 9l3 


923 


933 





V 












(31) 



We obtain the neutrino mass matrix as 

1 



mo 



Mr 



T 



Mr 







,2a 



/ e 

1 

V 1 

I 



{hri\){{9^i\r\{h^,} 



/ e 
1 
V 1 



Mn 



(32) 



As shown in Ref. ||2^, this mass matrix can naturally lead to a large v^-v^ mixing angle, 
which is suggested from the atmospheric neutrino oscillation ||2|. It is very crucial that 
the FN charges of the right-handed neutrinos are completely canceled out in the neutrino 
mass matrix in Eq. (|32|) and hence the masses of the neutrinos are determined only by 
the charges of the lepton doublets, (a, a, a + 1). This gives a mild mass hierarchy. 



m 



V2 



1 : 1 : 



C(l) : C(l) : C(10- 



(33) 



Now let us turn to our FN model. To change the above point, we suppose that the 
broken FN symmetry is not a U(1)fn but a discrete symmetry Z„ with n = 2d. Then, 
the mass matrix for the right-handed neutrino Mr changes into the following form; 



M, 



R 




9i2 e 

922 



c+d 



2c 



923 e 



b+c 



9l3(^ 
923 ^ 
933 e 



b+d 



b+c 
2b 



(34) 



Here we assume < b < c < d. Notice that the Majorana mass for A^^i is no longer 
suppressed by the power of e, which is a basic point to yield an extremely small neutrino 
mass m,yj. However, the structure of the neutrino mass matrix looks similar to the 
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original one; 




{hij} 



911 e-'''' 


9l2 


9l3 ^ 


-1 


/ e 







9l2 


922 


923 







1 




913 


923 


933 1 




lo 








(35) 



We see that one of the mass eigenvalue of this mass matrix strongly suppressed as 
~ g2(i+d)_ 'jj-^jg suppression is also understood directly by taking the determinant of the 
above mass matrix. For the mass hierarchy required from the successful AD leptogenesis, 
it is suitable to take d = 3 [Z^. 

To demonstrate our point, we randomly generate C(l) couplings hij and gij. Namely 
we calculate the mass matrix for neutrinos, taking the magnitudes of the couplings 
hij and gij to be in a range 0.5 - 1.5 and their phases to be - 27r. We also take 
e = 0.05-0.1 randomly.]^ Here, we have required the parameters r = bm^^^yjbm. 
« - )/K. - <), sin22^atm = 4|t/^3p (1 - \UM and tan^^sd = \Ue2lU, 



2 

atm 
2 



'el\ 



to 



be consistent with the parameter regions shown in Figs. ^ and ^. We have here required 
\Uez\ < 0.15 to satisfy CHOOZ limit P5[. Fig. R| shows the obtained mass of the lightest 



neutrino, m,^^. We can see that, an ultralight neutrino of mass rriiy^ ~ (0.1-3) x 10 ^ eV 
is naturally obtained. 



5 Discussion and conclusions 

In this paper we have performed a reanalysis on Affleck-Dine (AD) leptogenesis taking 
into account of all the relevant thermal effects. Then, we have pointed out that the 
baryon asymmetry is produced almost independently of the reheating temperature Tr 
and it is determined mainly by the mass of the lightest neutrino m,y^ in a wide range of 
the reheating temperature Tr ~ 10^-10^^ GeV. Notice that such reheating temperatures. 



~ 10^-10^^ GeV, are naturally realized in a large class of SUSY inflation models ^1 
p2| . This reheating-temperature independence is a very attractive feature of the AD 
leptogenesis. 

Furthermore, we have shown that the present baryon asymmetry predicts the mass 
for the lightest neutrino in a very narrow region, m,^^ ~ (0.3-1) x 10~^ eV. We have also 



proposed an explicit model based on a Froggatt-Nielsen mechanism ||T5| with a discrete 



^^A similar calculation was done in Ref. pq], where they adopted the U(1)fn model. 
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symmetry Zg, where such an ultrahght neutrino indicated from the AD leptogenesis is 
naturally obtained. 

Such a small mass of the lightest neutrino means that the mass parameter m^^y^ 
contributing to the Qv(3(3 decay is determined by the masses and mixings of two other 
neutrinos. Actually, we have shown that \my^y^\ can be predicted with high accuracy, 
by observable neutrino oscillation parameters such as 5m1^^, ^T^loh ^i^^ ^soi and Ue^- In 



particular, when the large angle MSW solution is the case and sin^ 9so\\j8ml^^ is relatively 
large, the value of \my^yj\ is predicted as \mu^uj\ ~ 10^^-10^^ eV, which may be testable 



at future 0i//9/3 decay experiments such as GENIUS |14 
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Figure 3: The parameter range we have taken for the atmospheric neutrino oscillation. 
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Figure 4: The parameter ranges we have taken for the large angle MSW, the small angle 
MSW and the LOW solutions to the solar neutrino problem. 
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Figure 5: The predicted value of \m,y^,yj for the large angle MSW solution. The solid 
lines represent the upper and lower values of \m,y^i,J for m^,^ ~ 0. The plots represent 
the values for the case when the m^,^ is allowed to be m^,^ < {1/ ^/2)mi^^. Here, we have 
required \Ue3\ < 0.15, to satisfy the CHOOZ limit. 
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sin^^soiT^mL [eV] 
Figure 6: Same as Fig. |, but for \Ue3\ < 0.10. 
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Figure 7: The predicted value of \mi,^„J for the small angle MSW solution. The solid 
lines represent the upper and lower values of Im^^^J for m,^-^ ~ 0. The plots represent 
the values for the case when the m^,^ is allowed to be m^-^ < (1/ \/2)m„2- 
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Figure 8: Same as Fig. but for the LOW solution. 
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